What a single snapshot reveals about the future and the past of turbulent flow 
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We develop an analytic formalism and derive new exact relations that express the short-time dis- 
persion of fluid particles via the single-time velocity correlation functions in homogeneous isotropic 
and incompressible turbulence. The formalism establishes a bridge between single-time Eulerian 
and long-time Lagrangian pictures of turbulent flows. In particular, we derive an exact formula for 
a short-term analog of the long-time Richardson law, and we identify a conservation law of turbulent 
dispersion which is true even in non- stationary turbulence. 
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Introduction. There are two alternative types 
of deseription in physies: one either uses a coordi- 
nate system fixed in space or follows particles. In 
fluid mechanics, the former is called Eulerian and 
the latter Lagrangian Respectively, our knowl- 
edge of turbulence roughly can be divided into two 
parts. On the Eulerian side, we have experimen- 
tal and numerical data on the single-time (snap- 
shot) velocity statistics: the moments of the ve- 
locity difference u measured at the distance R are 
(u") oc R^" . As for theory, the only analytic result is 
(uf) = —12eR/d{d+2) where e is the energy spectral 
flux, negative in 2d and positive in 3d, ui = {u-R)/R 
and d is space dimensionality 0, For d = 3, it 
is Kolmogorov's 4/5-law. On the Lagrangian side, 
the long-time growth of the inter-particle distance 
was inferred by Richardson from atmospheric data: 
{R^{t)) ~ \^\t^ , which was neither properly observed 
in a controlled laboratory experiment nor derived 
analytically. Here and below, all averages are over 
an ensemble of pairs whose separation initially was 
the same: .R(O) = Rq. Kolmogorov and Richardson 
laws, respectively Eulerian and Lagrangian, are ap- 
parently related by e and yet it is not known if the 
latter is a consequence of the former. 

Here we develop a formalism, which is a bridge 
between Eulerian and Lagrangian approaches. We 
derive a relation, which is a Lagrangian consequence 
of the Kolmogorov law and a short-time counterpart 
to the Richardson law. We express the Lagrangian 
time derivatives of the correlation functions of u and 
i? at t = in terms of the Eulerian moments of 
u. For short times, these determine the temporal 
evolution of the correlation functions by a Taylor 
expansion in t. Previously, analytic results on La- 
grangian evolution were derived either for spatially 
trivial smooth flows at long-time limit or for tempo- 



rally trivial short-correlated Kraichnan model. We 
derive some new results for these two cases, as well 
as for the general case of spatially non-smooth tem- 
porally correlated flow which corresponds to the in- 
ertial interval of turbulence. 

Of most interest is to flnd quantities conserved in 
time or at least having time derivatives all zero at 
i = 0; as we shall see below, the latter does not guar- 
antee the former. The only fact known before about 
two particles is that {[R{t)\~'^) does not change at 
i — > cx) in a spatially smooth flow (in the viscous 
but not inertial interval) One may think that 
the time dependence saturates only in the long-time 
limit when starting conditions are forgotten. How- 
ever, here we show that {R~'^) is an all-time genuine 
integral of motion for statistically isotropic spatially 
smooth random flow, even with a non-stationary 
statistics. To turn time derivatives into zero we shall 
use the simple mathematical fact: 



5[R) 



(1) 



Using that it is zero for i? 7^ 0, we find the general 
conserved quantity: (u''-™i?-''-"-"(u ■ ^)™). 

For the Kraichnan model with spatially non- 
smooth velocities we show that the moment, which 
generalizes (i?"''), has all time derivatives zero at 
i = 0, although it is not constant in time. For 
real turbulence in the inertial interval (spatially non- 
smooth yet temporally finite-correlated), we find 
that a direct analog is (^R^~^^~'^) . It has the first 
two derivatives equal to zero, but not the third one, 
showing what little traces of the conservation are 
left away from the limiting cases of velocities either 
spatially smooth or temporally white. This moment 
is nonetheless of much physical interest as both a 
short-time analog of the Richardson law, whose evo- 
lution is determined by the energy flux, and as a 
measure of the irreversibility of particle dispersion. 



statistics of the distances. Let us first cal- 
culate the time derivatives of correlation functions 
involving only R using the pdf 



V{R, t; Ro, 0) ^ (^S (^R - Ra - u{t')dt' 
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The integral is over Lagrangian trajectories and av- 
eraging here and everywhere is over u. To calculate 
time derivatives at t = of any mean, 



{F{R{t))) = j F{R)V{R,t-Ro,Q)dR , 



(3) 



we first differentiate the (5-function (and velocities), 
then set i = turning the (5-function into 5{R~ Ro): 
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where d/dt = dt + u-\I f( is implied. Here and below, 
unless stated otherwise, averages without an explicit 
specification of the time are taken at t = 0. 
When the velocity statistics is stationary, 

= (-^V.u,)V, , 
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Velocity time derivatives (acceleration etc) appear 
only starting from n ~ A. Therefore, in a general 
case the velocity snapshot contains enough infor- 
mation to describe evolution up to terms. For 
translation invariant incompressible flow, the first 
two derivatives can be simply obtained from the ve- 
locity single-time pdf f{u) assuming that the veloc- 
ity difference does not change with time: 

r{R,t;Ro,0) = \t\-^,f[{R-Ro)/t] . 

This ballistic regime was first considered in [HI, and 
it was recently suggested to be relevant for many 
physical situations [6|. 

Consider now general scalar correlation functions 
of the form {F {u{t) , R{t))) which may also include 
time derivatives d^u/df^. Let us label the particles 
by the indices 1, 2 and denote Eulerian spatial coor- 
dinates by r. Taking the Lagrangian time derivative. 
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at t = 0, one can set R ~ Rq in the term with 
the partial time derivative. It then turns into an 



Eulerian correlation function, where the average is 
done over all pairs in the flow separated by the same 
distance. Using incompressibility and homogeneity 
we rewrite the sum of spatial derivatives and obtain 



+ V,{u'F{u, Ro)). 



(5) 



We apply the formalism first to a spatially smooth 
random flow with isotropic but not necessarily sta- 
tionary statistics. In this case, = n i.e. the rel- 
ative velocity is proportional to the distance, u = 
d-{t)R{t), where (j{t) is a random traceless matrix 
with isotropic statistics. Such is the case, for ex- 
ample, in the viscous interval of turbulence, where 
the separation between the particles is much smaller 
than the viscous scale, for elastic turbulence and 
many other cases related to dynamical chaos. It has 
been shown by Zel'dovich et al. ^ that {[R{t)]-'^) 
does not change at t — >■ cx) in a steady flow, using the 
formalism of Furstenberg Q . Another proof under 
the same conditions is in [8| . Let us prove that all the 
time derivatives aX t = are zero for F{R) = R^'^. 
This is a consequence of (^d"'R~'^/dt'^) = when 
taken at R = Rq without additionally setting i = 
in ff{t). To show this we note that the first time 
derivative at R = Rq is zero by virtue of isotropy: 



(dR-ydt) = {<T,,(t))Rl,RlR„ 
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We can now use induction and express the n + 1 time 
derivative at R = Rq similarly to (l5|). 
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On the right side, the first term is zero by the induc- 
tion hypothesis, and the second term is zero because 
of (P) and (u*d"i?-''/(it") oc R^Rq"^ at R ^ Rq. An 
alternative derivation is given below by 



pUfTSj) for 

the stationary case. Our simple proof shows that 



([i?(t)]-'^)= V: 



(7) 



for all t even for a time-dependent statistics. The 
universality of this statistical conservation law fol- 
lows from the simple dynamical statement: in ev- 
ery realization of an incompressible spatially smooth 
flow, the integral of [R{t)]^'^ over the directions of 
the initial vector, Rq, is constant in time. That 
statement can actually be found by taking a closer 
look at the argument (not only result) of For 
isotropic fiows this integral can be interpreted as an 
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average over the angular degrees of freedom, which 
gives ([7]). Understanding that ([7]) holds for all times, 
even in the case of decaying turbulence, opens the 
door for its experimental verification, lacking so far. 

Non-smooth velocity. Now we ask if a simi- 
lar conservation law exists in a spatially non-smooth 
case. The only known results arc for the Kraich- 
nan model where the velocity is statistically sta- 
tionary, homogenous and delta-correlated in time: 
{u,{t)ujit')) = Dij{R)5{t - t'). The pdf satisfies 
the equation dtP = \7iDij\/jP = MV, see e.g. 

For an incompressible statistically isotropic flow, 
Dij{R) = R-^[{d+l--f)S,jR^-{2--f)RiRj], where 
7 is thus the measure of velocity non-smoothness, 
and time scales as R''. The role played by (-R^'*) 
in a smooth case is now assumed by (i?''^''), which 
follows from AdR"^^"^ = S{R). Indeed, this means 
that the first time derivative at any t is proportional 
to the probability of two particles initially separated 
by a distance Rq to come together at time t 0]: 



PiO,t;Ro,0), 



(8) 



Using the pdf from [8|, |9|, we derive: 
V{0,t-Ro,0) oc i?^-i|t|-''/^exp(-di?2/7'|t|) • (9) 

So, while (^R''^^'^ is not conserved in the Kraichnan 
model, all its time derivatives are zero at i = since 
V{Q,t; Ro,Q) has a substantial singularity. Time 
derivatives in the Kraichnan model differ from 
where we first set t = and only then average over 
velocities. The two procedures commute for finite 
correlated velocities but not for J-correlated ones. 

Let us see if there is any special moment F{R) = 
R^ in a real turbulence where velocities have finite 
temporal correlations. From now on we consider sta- 
tionary statistics. The first time derivative is zero 
for any h due to isotropy, then, assuming in addition 
incompressibility and translational invariance. 
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where ak = d — k + + b. For all moments but 
one, the second derivative is non-zero and short-time 
evolution is quadratic in time i.e. ballistic. 

The exceptional moment R^^^^^"^, which has both 
first and second derivatives zero, is a direct analog 



of i?'*'"'' in the Kraichnan model. The similarity be- 
tween these two cases is related to the fact that the 
operator {d^ /di^) — Vi{uiUj)\7j = A4 has the same 
structure as the operator A4 = WiDijVj, which de- 
termines the time evolution of the distance pdf in 
the Kraichnan model. Moreover, zero time deriva- 
tives at t = of (i?(t)^~^^~'') also originate from 
MR^^^'^~'^ = S{R). However, the third time deriva- 
tive (|10p of i?2-?2-d is non-zero, in distinction from 
the Kraichnan model. 

Similarity to the Kraichnan model may tempt one 
to conjecture that {d'^ R^'^^''^ /dt'^) cx P{0,t;Ro,0) 
for finite correlated velocities. To show that this 
is not the case, note that d]^V{R,t; Rq,0) at t = 
is proportional to S{R — Rq) that is zero for 
R / Rn. Thus, 7^(0, t; i?o, 0) has a substantial 
singularity at t = like in the Kraichnan model. 
Since (d^i?^~^^^'^/(it^) has a non-zero first deriva- 
tive, then it cannot be proportional to V{0, t; Rq, 0) 
in finite-correlated fiows. Indeed, it is quite general 
that a probability to cross a finite distance (from Rq 
to 0) has a substantial singularity at < = 0. On the 
other hand, non-analyticity of the moment evolution 
{R^{t)) is an artefact of the delta-correlated model 
and does not take place for a finite-correlated flow. 

Let us now describe the short-time evolution of 
([i?(f)]^~^2~'') using the energy flux relations in 2d 
and 3d combined with equation (fTO)l and the Kol- 
mogorov scaling ^3 = 1 and ^2 = 2/3. We get re- 
spectively for 2d and 3d: 
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(11) 
(12) 



In 3d, where 6 > and the energy flows to small 
scales, growth of this negative moment means that 
the main contribution comes from converging pairs. 
In 2d, where 6 < and the energy flows to large 
scales, diverging pairs dominate and the moment de- 
cays. These relations provide the Lagrangian conse- 
quence of the 4/5-law and the short-time counter- 
parts to the Richardson law. 

The moment (^R?~^'-^~'^^ thus provides an alterna- 
tive way of measuring 6, which is valid for short times 
and so should be accessible experimentally. Further- 
more, the initial dynamics of this special moment 
is irreversible in time, unlike other moments of the 
separation where irreversibility is hidden by the bal- 
listic evolution. At long times, we expect <^i?^~''^"'^) 
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to decay both in 2d and 3d, since diverging pairs 
should dominate the statistics. Indeed, if the pdf is 
a self-similar function of Rl^ /et^, which seems to be 
the case at least in 2d, we obtain {R'^/'^''^) oc ^2-3^/2 

We thus find that in the inertial interval of turbu- 
lence (spatially non-smooth velocity with finite tem- 
poral correlations), there are no moments whose first 
three time derivatives turn into zero; we conclude 
that the velocity snapshot completely determines the 
short-time evolution of the statistics of the distance 
between fluid particles. 

Velocity-distance correlation functions. For 
stationary statistics, dt{F{u,Ro)) ~ 0, so that ([S]) 
turns into a continuity equation: 

|(mt),i?(t)))^^^.^(.^n^^,i?o)) (13) 

The particular case is the Kolmogorov law expressed 
in a Lagrangian language: dt{u^) = (du^/dt) ~ 
'Vi{u^u'^) ~ — 4e [1, We can now generalize 

{du^R^/dt) = -4ei?g (1 + b/d) . (14) 

The initial evolution of the moments depending on 
the energy is controlled by e. The exception is as 
an effective energy density, {[u[t)Y[R[t)\~'^) , whose 
first derivative is zero and the short-time behavior 
is dominated by a time reversible contribution. 
More generally, {dufR''/dt) = cii?o"^(u°+^) and 
{d^ufR''/dt^) = C2{b + ac2 + a)i?^-2(M°+2)/(a + 1), 
where Ck = d — k + £,k+a + b. One can show that 
no a, b exist which turn both derivatives into zero in 
a non-smooth case. We thus cannot build statisti- 
cal integrals of motion compensating growth of u by 
that of R. Consider even more general form: 

G(m, R) = M"-'"i?^~"(u • RY" (15) 

with arbitrary a, 6, to. From our short-time perspec- 
tive, if {G{t)) doesn't have all time derivatives equal 
to zero at the initial time, it is surely not an integral 
of motion. Using stationarity and (jl3p . we express 
the condition for time derivatives of {G{t)) to vanish: 

The last equality would be satisfied by virtue of ^ 
if {u,d"'-^G/dt'^-^) = ConstRlR-^'^. If b is not a 



positive integer, that correlation function contains 
the term 

cx rI'"+^"^-RIj (17) 
for any n. Assuming that all terms in that function 
have the same scaling, we obtain the condition to 
satisfy (HH): 

6 - ?i + e„+a = -rf. (18) 

It is valid for all n only for the smooth case ^„ = n 
and b+a = — rf, which generalizes ([7]). We thus found 
a whole family of quantities whose time derivatives 
are all zero at t = 0. To claim conservation, we 
need analyticity in time, which holds when (jl5p is an 
analytic function of the velocity (which is an analytic 
function of time). We thus have an infinite family 
of conservation laws in a spatially smooth random 
flow. The entire family originates from the same 
dynamical statement as the special case 

We mention briefly the vorticity cascade in 2d, 
where ^„ = n, however it is straightforward to show 
that the logarithmic corrections to the scaling pre- 
vent the existence of the integrals in the form ()15p . 
Therefore, statistical integrals of motion for a pair of 
particles are found, so far, only for spatially smooth 
random flows (and Id compressible random flow 
where the inter-particle distance itself is the statisti- 
cal integral of motion [lH). For Burgers turbulence, 
one can also find the analog of (fT4|) with d = 1/3. We 
conclude that for non-smooth fiows, including the 
Kolmogorov scaling ^„ = n/i, statistical integrals of 
motion that exploit basic symmetries and ([1]) do not 
exist. Either the integrals appear at the long-time 
limit or finding them requires deeper insight into the 
interplay between dynamics and geometry. 
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